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We point out that a light scalar field fiuctuating around a symmetry-enhaced point can generate 
large non-Gaussianity in density fluctuations. We name such a particle as an ^^ungaussiton" , a 
scalar field dominantly produced by the quantum fluctuations, generating sizable non-Gaussianity 
in the density fluctuations. We derive a consistency relation between the bispectrum and the 



c/3 . 
I ^i ' trispectrum, tnl ~ 10'^ "/nl ' '^bich can be extended to arbitrary high order correlation functions. 

^ ■ If such a relation is confirmed by future observations, it will strongly support this mechanism. 
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I. INTRODUCTION 

The concept of symmetry has been a guiding principle in modern physics. The structure 
of the standard model (SM) is dictated by the SM gauge symmetries, SU{3)c x SU{2)l x 
U{1)y- In the string theory, gauge symmetries as well as discrete global symmetries are 
ubiquitous [l|]. 

One of the plausible candidates for the theory beyond SM is supersymmetry (SUSY), 
which relates a particle to a superparticle with different statistics. In particular, a super- 
partner of a fermion field with spin 1/2 is a scalar field. Therefore there are perhaps many 
fundamental (or composite) scalar fields in nature, and most of them may be charged 
under some symmetries. 

A scalar field that is singlet under any symmetries does not possess a special point 
in its field space, and so, the minimum of the potential may vary in the evolution of 
the universe. The change of the minimum generically induces coherent oscillations of 
the scalar field. If the scalar field has only interactions suppressed by the Planck scale, 
Mp(= 2.4 X lO^^GeV), it might induce notorious cosmological modu 
has recently turned out to be more acute than previously thought P, \M, 



i^problem [2[, which 
a]. That is to 



say, the modulus decay generically produces too many gravitinos, which either spoil the 
success of the big bang nucleosynthesis or exceed the observed dark matter abundance. 
The moduli problem and the moduli-induced gravitino problem can be ameliorated if the 
initial displacement is somehow suppressed. This, however, is difficult to achieve unless 
the modulus field is charged under some symmetries. 

On the other hand, a scalar field charged under some symmetries has a special point 
in its field space ^, where the symmetries are enhanced and the scalar potential has its 
extremum. If the symmetry-enhanced point is the (local) minimum of the potential, it is 



quite likely that the sea 



continues to sit there 



ar field sits at the point from the beginning of the universe and 



y, y, y, y 



10|. Such a scalar field apparently does not affect the 



One exception is an axion with a shift symmetry, because the axion docs not have a symmetry-enhanced 
point. We will make a comment on this case later. 



subsequent evolution of the universe, since its energy density vanishes, at least classically. 
At the quantum level, however, this may not hold as we describe below. 

The inflation is now strongly suggested by the recent observations of the cosmic mi- 
crowave background (CMB). During inflation, the inflaton acquires quantum fluctuations, 
which become the seed for the present structures of the universe. Interestingly, any scalar 
fields of masses smaller than the expansion rate during inflation, if/, similarly acquire 



quantum fluctuations of 0{Hi) |lll ]. 



We point out that such light scalar flelds fluctuating around the symmetry-enhanced 
point necessarily induce non-Gaussianity at some level in the observed CMB power spec- 
trum ^. In other words, the scalar flelds whose energy classically vanishes must have 
left their traces in the CMB; those scalar flelds are produced by the purely quantum 
fluctuations, making the fluctuations in their energy density deviate from the Gaussian 
distribution. We call such a scalar particle an ^^ungaussiton" , because their energy density 



exhibits strong non-Gaussianity 1^, ll3|. Our deflnition of the ungaussiton is an elemen- 
tary or composite scalar field dominantly produced by the quantum fluctuations, generating 
sizable non-Gaussianity in the density fluctuations. Note that the deflnition itself is not 
related to the symmetry. The presence of the symmetry is one of the ways to naturally 
make a scalar fleld to be an ungaussiton. 

If there are indeed many scalar flelds with symmetries in nature, some of them may 
be light during inflation. Then they can be ungaussitons. In this way we expect that 
the presence of non-Gaussianity might be quite generic in our vacuum. Interestingly, it 
has been recently reported that large non-Gaussianity was detected by the analysis on 
the WMAP Syr data [15|. The latest WMAP Syr data seem to have the same tendency. 



although the vanishing non-Gaussianity is allowed within 95% C.L. [16|]. Those hints on 
the non-Gaussianity, or its future detection, may be originated from such particles, i.e., 
ungaussitons. 



^ As will be discussed later, the symmetries are not necessarily exact. Our arguments remain intact 
as long as the deviation from the potential minimum is small (or comparable) relative to the Hubble 
parameter during inflation. 



In this paper we study the ungaussiton mechanism and write down a condition for 
the non-Gaussianity to become large enough to be observed by the ongoing and planned 
observation such as the WMAP and Planck satellites. Moreover, we derive a consistency 
relation between the bispectrum and the trispectrum, which, if confirmed, will be a smok- 
ing gun of the ungaussiton mechanism. Importantly, the bispectrum and the trispectrum 
predicted in the ungaussiton scenario exhibit a specific dependence on the scale; both get 
enhanced as one goes to smaller scales. Such a feature on the scale dependence as well 
as the consistency relation between the bispectrum and the trispectrum will help us to 
distinguish our scenario from the others. 

Let us here comment on how the ungaussiton can remain around the minimum of the 
potential during inflation. As is well known, if inflation lasts sufficiently long, a light scalar 
field (an ungaussiton in our case) relaxes towards the Bunch-Davies (BD) distribution |l4| . 
The BD distribution is such that the root mean square (RMS) of the field fiuctuations on 
large scales is given by 0{H]/m), where m is the effective mass of the ungaussiton during 
infiation. If m is much smaller than Hj, the BD distribution tells us that the ungaussiton 
takes a value of 0{H]/m) ^ Hj in most regions, while there is a small portion of the 
whole universe where the field value is of 0{Hi). If our observable universe happens to 
be inside such small region, the spatial average of the ungaussiton over the observable 
universe can be very close to the origin. This may sound a fine-tuning, but it may not. 
There may be anthropic arguments to select such regions, although there is no justification 
for such reasoning at this moment. There are other ways to keep the ungaussiton near the 
origin. One possibility is that the mass m is not much smaller than Hj. Then the RMS 
values of the ungaussiton is not that large compared to Hj, which makes it easy to find 
such regions that the field value is of 0{Hi) without fine-tunings. Another is that the last 
infiation was not long enough. For instance, the e-foldings of 50 ~ 60 is too short for the 
field to reach the BD distribution. The ungaussiton then stays near the minimum of the 
potential throughout the infiation, if the ungaussiton dynamically relaxed to the origin 
before the inflation. We will come back to this issue in Sec. 



II. NON-GAUSSIANITY FROM AN UNGAUSSITON 

Let us consider a scalar field, a, which is charged under some symmetries ^. We set 
the origin of a to be the symmetry-enhanced point, which is assumed to be the (local) 
minimum of the potential. If it is a meta-stable vacuum, the life time of the vacuum is 
assumed to be much longer than the present age of the universe. The potential around 
the origin can be obtained by expanding the scalar potential in terms of a. We assume 
that the potential is well approximated by the quadratic potential *^, and that it is valid 
at least up to cr = 0{Hi): 

U{a) = \mla\ (1) 

where mo- denotes the mass of a °. The symmetry of U{a) is the Z2 symmetry in this case. 
If rrifj is larger than Hi, the fluctuations of a get suppressed, while a acquires fluctuations 
if mcr < Hi. In the following we focus on the case of m^ < Hi. 

We first give a rough sketch of how non-Gaussianity is generated in the primordial 
curvature fluctuations. During inflation, the scalar field a likely sits at the origin, a = 0, 
based on the symmetry argument. For the moment we assume this is the case, and we 
will discuss this issue in detail in Sec. IIVI Since the mass is lighter than the Hubble 
parameter, a acquires quantum fluctuations, 5a, around the origin. The magnitude of 
5(7 is roughly equal to Hi. After inflation, a starts to oscillate around the origin when 
the Hubble parameter becomes comparable to mg-, with an amplitude taking a different 
value in each local patch of the universe. The energy density of a is proportional to 
the amplitude squared, {5aY ^- Taking 5(j as a Gaussian variable, the resultant energy 



density exhibits strong deviation from the Gaussian distribution 12, 



ITj, hence the name 



'^ It is easy to extend our results to a case of the multiple scalars. 

'^ One can also consider a case that a quartic potential dominates over the quadratic one. The resultant 

non-Gaussianity tends to be smaller in this case. 
** As mentioned before, the effective mass during inflation may be different from that in the low energy. 

Our arguments remains intact as long as the effective mass is smaller than the Hubble parameter. 
Precisely speaking, one has to take account of the spatially averaged value of cr, which does not vanish 
in the observable finite universe. Then the dependence on 5a becomes slightly modified. Nevertheless 
this intuitive explanation gives a qualitatively valid picture. 



an ungaussiton. For simplicity, we assume that the ungaussiton a decays into radiation in 
the visible sector ^. Then, the fluctuations in the energy density of the ungaussiton turn 
into small corrections to the adiabatic fluctuations arising from the inflaton. Although 
the corrections might be small in the amplitude of the curvature perturbations, they can 
significantly change higher order correlation funcitions such as the three-point function 
which characterizes the strength of non-Gaussianty. 

Let us now estimate non-Gaussianity produced by the ungaussiton. First we define the 
fluctuations of a, taking account of the fact that our observable universe is finite. Assume 
that an ensemble average of a vanishes. Since our observable universe is finite, the spatial 
average of a over the comoving volume corresponding to the present Hubble horizon does 
not vanish. In order to calculate the spatial average of a, we need to define the time 
slicing. For later convenience, we take a fiat time slicing at t = t* when the smallest scale 
of interest leaves the horizon during inflation. Let us denote the spatial average of a on 
the flat slicing by o and define the fluctuation 5cr=i,(x)as bo^,[x) = cr{t^,,x) — a. If the mass 
of a is much lighter than Hj, those fluctuations defined on the slicing can be evaluated 
when each scale leaves the horizon. In the following we assume that the inflation lasts 
long enough that a acquires an almost scale-invariant Gaussian fluctuations which extend 
beyond the present Hubble horizon. To put it more precisely, we define the fluctuations in 
a box of a size L, which is (at least) several times larger than the present Hubble horizon. 
In fact, the fluctuations of the scales beyond the current horizon size contribute to the 
non- vanishing spatial average of a. 
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According to the 6N formalism [19|, |20|, |2l|, |2^ , on sufficiently large scales, the curvature 
perturbation ( on the uniform energy density hypersurface at t = t/ is equal to the 
perturbations in the number of e-foldings between the uniform density slicing at t = t/ 



s The ungaussiton can generate sizable non-Gaussianity even if it does not decay. In this case, the 
non-Gaussianity arises from the isocurvature fluctuations, which leave the characteristic signatures in 



the CMB anisotropy [18j . Indeed, the Z2 symmetry of the ungaussiton may naturally account for the 



stability to become dark matter. 



and the initial flat slicing at t = t*: 

C{tf,x) = N{tf,U,x)-N{tf,Q (2) 

with 

N{tf,t^,x) = / H{t,x)dt, 
Jt, 

N{tf,Q = r H{t)dt, (3) 

Ju 

where N(tf,t^,x) is the number of e-foldings of expansion between the two slicings, and 

N is that of the background universe (or more precisely, the spatially averaged one). 

Similarly, H{t,x) denotes the local expansion rate, while H{t) denotes the background 

one. 

We assume that the ungaussiton a decays after the reheating ^, and set tj at the time 

well after the decay of a. As mentioned, we take t^, at the time when the smallest scale 

of interest crossed the Hubble horizon during inflation. The curvature perturbation ( is 

conserved after t = tf, that is, ( becomes independent of t/, since the adiabatic pressure 

condition is satisfied. Then we can regard the number of the e-foldings as the function of 

0^ = (f)[t^,x) and cr* = a{t^:,x), and expand ( in terms of the fluctuations (50^, and 6a^:, 

which are the perturbations of 0* and a*, as 

C ^ N^Scf), + NJa, + -N^Jal + -N^^Ja't + —N^aaJ(Tt + ■■■ , (4) 

where the indices denote the partial derivative with respect to the variables, i.e., A^^ = 
dN/d(f), No- = dN/da, etc.. Those coefficients depend on the background evolution of the 
scalar fields </> and a from t = t^, until t = tf. In particular, we use the spatially averaged 
value, 0" = 0", as the initial condition of the background evolution. Here and in what follows 
we omit the dependence on the comoving coordinate x unless otherwise stated. We assume 
that the first term on the right hand side (RHS) of Eq. (jl]), coming from the infiaton 



'^ This is just for the sake of simphcity. If it is the infiaton that decays later, we should take tf after the 
reheating. However, the non-Gaussianity tends to be suppressed in this case, since the energy density 
of the ungaussiton becomes small. 

7 



fluctuations, dominates the two-point function of (. As for the relevant higher order 
correlation functions, we assume that the leading contributions come from the ungaussiton 
fluctuations, and we have neglected second and higher order terms in 50* in Eq. (j4]). It 
is known that the non-Gaussianity of ( coming from the intrinsic non-Gaussianity of 50* 



and 5cr* is far below the observational sensitivity such as Planck 23|, UM, l25|, l26| . Hence 
we can treat both 50* and 5o"* as the Gaussian variables '. 

Before going further, it will be worth mentioning how the symmetry comes into the 
game. The expression @ itself is rather generic, but, what is peculiar to the ungaussiton 
mechanism is the relative size of the coefficients. To see this, let us concentrate on A'^,., 
which depends on the background evolution of a. As mentioned above, we take a = a as 
the initial condition. Let us expand N^r around another initial condition, a = 0. 

Na\^ = N^\o+N^a\oa+-N^aa\o^^ + ---, (5) 

where we have explicitly shown the initial values of a. What the Z2 symmetry tells us 
is that the potential, and therefore the background evolution, is even around the origin. 
That is to say, the functions N^ and Ao-o-o- should vanish when a = is used as the initial 
condition. That is, 

A^.lo = 0, Ar,..|o = 0. (6) 

Therefore the expansion of N„\- around cr = actually starts with a term linear in a. 
Noting that a is actually comparable to the fluctuations 5cr*, one can see that the term 
linear in 5(T* in (jlj) is comparable to the term quadratic in 5(T*. Intuitively speaking, since 
a stays at the origin based on the symmetry arguments, if one expands a function of a, 
f{(T), around the origin, the leading term should start with cr^, as long as f{a) satisfles 
/(cr) = /(— 0") in the interested range of a. When one takes account of the flniteness of the 
current horizon scale, this feature appears in the suppression of N^\- as shown above. Let 
us emphasize again that the suppression is due to the presence of the symmetry, and that 



' Although we have defined Sa as the deviation from a, one can still regard Sa^ as the Gaussian fluctu- 
ations. This is because the superhorizon fluctuations contribute to a, while Sa represents fluctuations 
inside the horizon. 



this makes the higher order term quadratic in Sa,, very important. For the same reason, 
the cubic term of 5a^. is comparable to the quartic one, and therefore we keep the quartic 
term of 5(T^ for the moment. 

The power spectrum P^, bispectrum B,^ and trispectrum T^ of the curvature perturba- 
tions are defined as 



\3 



{CtSkM.)^ = W^ih + h + h) B^ih, h, k 



3 J 



\3 



{CkMUSk,)c = Wm + k2 + h + h) Tc(/ci, h, /C3, h), (9) 

where C^ are Fourier components of (, i.e., C^ = j d?xe~^^'^C,{x), ki = \ki\ for i = 1 ■ ■ -4, 
and the subscript c means that we take the connected part of the corresponding correlator. 
It is useful to define V(^ = A;^/(27r^) P<^. Let us also define the power spectra of 6(f) and 6cr, 

(%,%,) = i^nfm + k2)P4ki), (10) 

i^^kMk,) = i2n)'5ih + h)PAki), (11) 

k^ 

^0W = ^^-/-H- (^2) 



Neglecting the tilt of the power spectra, we have 



H,-' 



n-P,«(^j. (13) 

As mentioned above, the infiaton dominates the two-point correlator of (, i.e., 

n ^ KP^- (14) 

As for the bispectrum, the dominant contributions come from the ungaussiton, which, up 
to sixth order in 6a^, gives 

B^iku k2, h) ~ -^^ {P^{h)P^{k2) + 2 perms.) 

+ ^ / (0 ^C(^)n(l^i + ^)Pd\k2 - q\) 

+ ^"^^fr^ (/ ^ Pd<i)Pd\k2 - q\)Pdki) + 5 perms.) , (15) 



where the third term was not taken into account in Refs. 17|, |28|, |29|, ISOj. The first term 
is quartic order in 5a ^, while the remaining terms containing the momentum integral are 
sixth order in 5a ^. 

In a similar way, the trispectrum up to eighth order in 5a^: can be written as 

T^(fci, 4, ^3, 4) ^ ^^^ (Pc(A:i)Pc(^2)n(^i3) + 11 perms.) 
+ ^^ mh)P^{k,)P^{k,) + 3 perms.) 
+ ^ (/ T^ Pd^)PSki - q\)Pd\k2 + q\)Pd\ki + 4 - g1) + 2 perms.) 

+ ^-y— (I ^ Pc(?)n(l^2 + gDPcd^i + ^4 - g1)Pc(^i) + n perms. 

+ ^ ;^^8 '^'^'^ (y (^ n(?)n(l^3 + g1)Pc(^i)n(^i2) + 11 perms. J 

+ ~^2N^ [j W? ^^^^^^^^'^^ + ^' + g1)Pc(^^i)n(^2) + 11 perms. j 

+ ^-^-^— (I ^ Pc(?)n(l^3 - g1)Pc(^i)n(^2) + 11 perms.) , (16) 

where kij = \ki + kj\. We have here dropped small logarithmic corrections proportional to 
og (L/i?), where R denotes a smoothing scale. The first three terms agree with those in 



17 



29 



30|. 



It is conventional to express Br and Tc in terms of the non-linearity parameters /nl, Tnl 



and gjsiL defined by |27| 



Bciki, k2, ks) = ^hL{Pdki)Pdk2) + 2 perms.), (17) 

T(;{ki, k2, h, h) = ^nl(^c(^13)^c(^3)^c(^4) + H perms.) 

54 
+—9MPdki)Pdk2)Pdh) + 3 perms.). (18) 

If we include the terms containing the momentum integral in i?^ and T^, the parametriza- 
tion of Eqs. (fT7I) and (fT8|) does not correctly represent the bispectrum and the trispectrum. 
However, if we consider some limiting configurations of the wavenumber vectors, so-called 
"squeezed" configurations, we can approximately express B^ and T^ in terms of /nl, tnl 
and fifNL- 

10 



Let us here make another approximation which drastically simplifies our arguments. In 
the ungaussiton scenario, it turns out that those terms containing A^^^-o-o- or Ncraaa in T(^ 
are suppressed compared to the terms containing only A^^- and N„„, by the ratio of the 
ungaussiton energy density to the total one at the time of its decay. The ratio must be 
smaller than ~ 10~^ not to exceed the primordial fiuctuations produced from the infiaton. 
Hence if some quantity of interest starts from terms which contain only A^^- and/or N„„, 
we will drop higher order terms including A^o-o-o- or N„„„„. 

In the case of the bispectrum, the squeezed configuration is such that one of {ki} with 
z = 1, 2, 3 is much smaller than the length of the other two wavevectors; e.g., ki <C /c2, k^. 
For such configuration, /nl can be written as 



17 



28 



29 



3Q| 



^/nl ^ ^ (A^.'A^.. + NlP, logihL)) , (19) 

where kf, = min{fcj} {i = 1, 2, 3) and L is the size of a box in which the perturbations are 
defined. Since we are interested in the perturbation in the observable universe, L should 
be taken to be 0{1/Hq), where Hq is the present value of the Hubble parameter. 

As for the trispectrum, the situation becomes a littile bit more complicated. If the 
four wavenumber vectors are such that minj/cj, \kj + ki\} is much smaller than the other 
elements of {/cj, I fcj + A;^ I}, where {i,j,i) = {1, 2, 3,4}, tnl can be expressed by IT], |29|, |30| 

^NL ^ ^ {NXa + KaT. log(fc,L)) , (20) 

— * — # 

where kt = in.m{ki, \kj + ki\}. Meanwhile, if the configuration of the four wavenumber 

vectors is such that two of them is much smaller than the other two, (^nl is given by 

54 1 

— ^NL ^ ^ {N'^N^^^ + mlN^^N^^^^VAogikgiL) + m^Nl^N„,^VAog{kg2L)) , (21) 

where k^^ = minj/cj/c^fc^} and kg2 = TLn.in{kfkj}. Note that (^nl starts from the terms 
linear in A^o-o-o- or N^^aa- Hence the trispectrum is dominated by tnl terms if the leading 
non-Gaussianity is generated by the ungaussiton. Here and in what follows we neglect g'NL 
terms and the tilt of the power spectra, for simplicity. 

The non-linearity parameters depend on the thermal history of the universe after in- 
fiation, since the coefficients such as N„ and N^a depend on the evolution of the infiaton 

11 



and the ungaussiton. We consider a case that the ungaussiton decays after the reheating, 
since we are interested in the case that the non-Gaussianity becomes large. That is to say, 
we focus on the case of 

r. < r,^, (22) 

where ^a{(t>) is the decay rate of the ungaussiton (inflaton). There are still two cases to be 
considered, depending on whether the ungaussiton starts oscillations before or after the 
reheating. 

If the ungaussiton starts its oscillations after the inflaton decay, or equivalently, if 
rricr < r^, a relevant part of the e- folding number A^ that depends on a is given by 31 1 



Nitf,Q D ^^a^p + Oip 



2\ 



(23) 



where a = 2y^2/'nV{b/A) ^ 1.45 is a numerical constant and p = {a j Mpf \fniaJ^cT- The 
parameter p roughly represents the ratio of the ungaussiton energy density to the total 
one at the time when the ungaussiton decays. If the ungaussiton dominates the universe 
before its decay, the resultant curvature perturbations would become highly non-Gaussian, 
which is inconsistent with the current observations. Hence we require p <^ 1. Note that we 
have assumed that the curvature perturbation dominantly comes from the inflaton, while 
the ungaussiton gives only small correction to it. This actually gives a severer bound on 
p, i.e., p < O(10~^). Substituting ^ into ^ and ([20]), we obtain 

where e is one of the slow-roll parameters defined by e = A/p/2 jV^^/Vp iy is the scalar 
potential of the inflaton) , and it is related to A'^^ and "P^ as 

The WMAP normalization gives "P^ ^ 2 x 10"^ [id]. 

Let us here discuss the magnitudes of the terms in the curly bracket of (!24l) and (l25l) . 
Since we are interested in a scalar fluctuating around the origin, the first term should be 

12 



0(1). There are several ways to realize such situation as stated in Sec. HI In order to 
estimate the magnitude of the second term, we need to specify values of k and L. What 
values of L should be taken depends on the spatial size of the observational data. For 
instance, suppose that we have observational data on a variable X in a small patch of the 
sky. Then we should take L to be the size of the patch. We can define a homogeneous 
part of X by taking the spatial average over the region. The fluctuations is obtained 
by subtracting the homogeneous mode from X. Those fluctuations with the wavelengths 
larger than L cannot be ditinguished with the homogeneous part, and therefore they are 
regarded as the homogeneous mode. In our case, since we are interested in the fluctuations 
of the scales between k ^ Hq and k ~ lO^ifo, we should take L ~ 0{Hq^), and the second 
term is 0(1) for interested ranges of the scales. Therefore, both terms in the curly brackets 
are 0(1). 

The fact that both terms in the curly brackets are 0{1) means that if the ungaussiton 
is responsible for the non-Gaussianity of the curvature perturbations, the non-linearity 
parameters /nl and tnl have the logarithmic dependence on the spatial scales K In 
particular, they increase logarithmically as one goes to small scales. 

When does the non-Gaussianity become large? The condition, /nl > fm^, is met if 

(p \ —3/2 / rmin \ 

where we have set the terms in the curly brackets on RHS of (!24|) and (125|) to be unity. 
Equivalenty, we can express the same condition in terms of the Hubble parameter during 
inflation instead of the slow-roll parameter e; 

/n \1/V T \1/V ^ \ -1/2 /fmm\l/6 

U2 r^ \r I y*^ 1 I (^ \ I °^^ \ I •'NL 



where T^ and Tqsc are the temperatures of the radiation when the ungaussiton decays and 
when the ungaussiton starts to oscillate, respectively. Here g^^i and g^^osc are the effective 
number of light degrees freedom at T = T^ and Tqsc, respectively. 



J Note that we have neglected the scale dependence oiV(^. However, it can be determined by observations, 
and so, in principle one can extract the scale dependence of the non-linearity parameters. 

13 



Let us next consider the other case that the ungaussiton starts its oscillations before 
the inflaton decay. In this case, we have only to replace rricr with F^ in Eqs. fl27|) and fl28l) . 
employing the sudden decay approximation on the inflaton decay. Hence the condition, 
/nl > /n™) can be written as 






min 



e'if >7xl0^"(^), (29) 



or equivalently. 



(\ 1/8 / rp \ 1/2 / T-i \ -1/2 / j-min\ 1/6 

S) (r^) iJkv) (if) ■ (») 

where Tr is the reheating temperature of the inflaton, and g^,R the effective number of 
light degree of freedom at T = Tr. 

There is the upper bound on the Hubble parameter during inflation as Hi < 10^'^GeV 



from the WMAP result IGj. Meanwhile, both (l28l) and (l30l) give the lower bounds on 
Hj. Notice that both the conditions (128|) and (130|) require the decay temperature be much 
lower than the reheating temperature, i.e., the life time of the ungaussiton must be rather 
long. The longevity is necessary for the ungaussiton to give non- negligible contribution to 
the energy density, and therefore to the non-Gaussianity. It suggests that the ungaussiton 
must couple only weakly to the visible sector. Suppose that a couples to the visible 
sector strongly. Then a can be immediately thermalized when the universe is reheated 
by the inflaton decay, which would signiflcantly suppress the non-Gaussianity. Thus the 
ungaussiton must be indeed weakly coupled to the visible sector. 

In Fig.lll we show the contours of the non-linearity parameters /nl and tnl as functions 
of Td/Tosc (or Td/TR) and Hi. In the flgure, we set ^f^^osc = g*d = 9*r = 100 and a = Hi/27i, 
and drop the logarithmic dependence for simplicity. We see that the region of large 
Hi and small T^/Tosc (or Td/Tn) is already excluded by the observational constraint. 



|/nl| ^ 100 [l6|. We also flnd that the slopes of the contours do not depend on the non- 



linearity parameters: Hi oc a/T^^/Tosc oi^jTd/Tn for all the non-linearity parameters. This 
is because both /nl and tnl depend only on the combination e^fm^Jvl (or e-^/T^/TV), 
as can be seen from Eqs. (l24l) and (l25l) . 
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FIG. 1: Contours of the non-linearity parameters /nl and tnl; (/nLi '^nl) = (100, 5 x 10^) (solid 
line) and (1, 10^) (dotted line). Shaded region is already excluded by the WMAP 5yr data. Here 
we set 5*osc = g*d = g*R = 100 and a = Hj/Itt, and drop the logarithmic dependence. 



The fact that all the non-linearity parameters depend only on one combination of the 
model parameters means that we have an universal relation between /nl and Tnl, which 
is independent of the model parameters. Such a consistency relation is given by 



T^L = CV,' (-/nl 



C 



1 X 10=^ C /^/' 



NL ! 



(31) 



where we have used Pj; ~ 2 x 10 ^ 16|. Here C is a numerical coefficient of order unity 



15 



defined by 



—^—7 + log hL 
C = ^ (32) 



,4/3- 

^ + log hL 



As a remarkable fact, when /nl > 1, the trispectrum of the curvature perturbations be- 
comes inevitably very large tnl ^ 1- Although the current bound on tnl is very weak, we 
will have much stronger bound in the near future: e.g., ItnlI ^ 560 for Planck 32|. Hence 
if the non-Gaussianity /nl > 1 is generated by the ungaussiton mechanism, we should 
also see the strong non-Gaussianity through the trispectrum, which will be quite useful 



to distinguish the ungaussiton scenario from the others that 



relations between the bispectrum and the trispectrum 31 



Dredict different consistency 



33 



m- 



Let us here show that the consistency relation (|3T1) can be understood in an intuitive 
way. The bispecturm is roughly given by {C,'^) ~ A'^^^ ((5cr)®), while the trispectrum is 
(C^) ~ A^^ (((5cr)*). Therefore, as long as the quadratic term (5cr)^ is the dominant source 
of the non-Gaussianity, we naively expect (C^) ~ ((C'^))^^^) which correctly reproduces 
Tnl ~ '^c~^ /nl ~ 10'^ /nl • ^^ '^^'^ easily extend this argument to the ra-th correlator. 
The leading n-th non-linearity parameter defined in a similar fashion, f^^ , is roughly equal 
to Vc; f^i^ ~ 10'^*-"~'^''/^L • 0^6 can see that this simple argument breaks down if Ao- 

is not suppressed as in the curvaton scenario, since the bispectrum would be dominated by 
{(^) ~ NlNfjfj ((i^cr)^)- In that case there are different consistency relations [3l|. Thus it 
is very important to test the consistency relations by observations, in order to distinguish 
different scenarios to produce non-Gaussianity. 

As mentioned before, there must be hierarchy between Tr (or T^sc) and T^. Since the 
ungaussiton must decay before the big bang nucleosynthesis (BBN), Td must be higher 



than 10 MeV 



351]. Then, the reheating temperature is bounded below: e.g., Tr > 10^ GeV 



for Hi = lO^^GeV, and Tr > lO^GeV for Hi = 10^^ GeV, in order to have /nl > 1- 
Such high reheating temperature as Tr > lO^GeV may lead to the overproduction of the 
gravitinos |36|, |37| , when the local SUSY is assumed. One way to ameliorate the tension is 
to consider the heavy or ultralight gravitino mass: m3/2 > 0(10) TeV or m^/^ ^ 0(10) eV. 



Another is to introduce a late-time entropy production after the ungaussiton decays [38 
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Note that the density fluctuations after the ungaussiton decays are adiabatic, although 
they contain sizable non-Gaussianity. Once it becomes adiabatic, the non-Gaussianity will 
be inherent in the subsequent evolution of the universe, whatever it is like. In particular, 
the non-Gaussianity does not change even if there is late-time entropy production after 
the ungaussiton decays, unless isocurvature perturbations exist in the flelds that induce 
entropy production. 

Let us here make a comment on the isocurvature perturbations. If the baryon number 
or the dark matter abundance is flxed before the undaussiton decays, the non-Gaussian 
isocurvature perturbations are generated. The amplitude of the isocurvature perturbations 
are smaller than the current bounds, for the non-Gaussianity with /nl < 100. Although 
we do not consider its effects throughout this paper, the isocurvature perturbations, if 
observed, may provide another clue on the ungaussiton mechanism. 

Before closing this section, we note that the symmetries does not have to be exact to 
realize the ungaussiton scenario. The explicit breaking of the symmetries will generically 
makes the potential minimum deviate from the origin, and a may have a large value. 
However, if the deviation is the order of if/, this is indistinguishable from the non- vanishing 
a due to the fluctuations on large scales beyond the size of the observable universe. Also, 
as in the case of an axion fleld with a shift symmetry: a — > a -|- a, the enhanced-symmetry 
point may not exist. Even in this case, if the initial position of the scalar happens to be 
near the minimum, it can be an ungaussiton if it is light during inflation. Therefore our 
discussions above remain intact in these cases. 

III. MODELS 

Let us here give some models to realize our idea. We will give detailed analysis on 



each model in the coming paper 40|. One example is the modulus fleld, which appears in 
the supergravity/superstring theory and it has a very flat potential. If the position of the 
modulus field during inflation is deviated from the minimum in the low energy, it leads to 
the cosmological moduli problem, since the large energy density stored in the modulus fleld 
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can easily modify and spoil the standard evolution of the universe. One of the solutions 
to the moduli problem is to presume that the moduli fields having symmetry-enhanced 
points in their field space stay at the points during and after inflation [y, |7|, |8|, y, [lO| , which 
would greatly suppress the cosmological abundance of the moduli. Such moduli fields 
become ungaussitons, if they are light during inflation ^. Since the moduli must decay 
before BBN starts, the modulus mass must be larger than O(IO^) TeV. Assuming that the 
modulus has couplings with the standard-model particles suppressed by the Planck scale, 
the decay temperature is given by around O(IOO) MeV for the modulus mass of 10^ TeV. 
From Fig. [1], one can see that large non-Gaussianity is generated for e.g. Hj = 10^^ GeV 
and Tfj = 10^'^ GeV. In the case of the moduli, the needed hierarchy between Tr and T^ is 
realized due to the weakness of the modulus couplings. 

Another is a right-handed sneutrino, A^, a superpartner of the right-handed neutrino, 
A^. The right-handed sneutrino is odd under the -R-parity, which determines the origin 
of A^. If the mass of A^ is light during inflation, it can be a ungaussiton. If the neutrino 
Yukawa coupling y is small enough and if the mass of A^ is relatively light, the decay of 
A^ can be delayed. Also one can neglect the thermal production of A^ for sufficiently small 
Yukawa couplings. Therefore, it is relatively easy to satisfy the relation ( !3Ql) if the model 
is close to the case of the Dirac neutrino. For instance, the decay temperature is roughly 
estimated by 



Td ~ 0.1 y^m^Mp, (33) 

where m^ is the mass of A^. Therefore, large non-Gaussianity is generated for e.g. 
y = 10-^ mjv = 1 TeV, Hi = 10^^ GeV and Tr = W GeV. Note that the gravitationally 
produced right-handed sneutrino cannot explain the baryon asymmetry through leptoge- 



nesis |42[, since the baryon asymmetry would have too large isocurvature perturbations. 
The last example is a flat direction of the supersymmetric standard model (SSM). 
There are many flat directions in SSM, and they are parametrized by a gauge-invariant 
monomial such as udd and eLL. The flat directions are composed of squarks, sleptons and 



^ The abundance of the such modulus field was investigated in Rcf. 4l| to study how the moduli problem 
is relaxed by the mechanism. 
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Higgs, and therefore, the origins of the flat directions are determined by the SM gauge 
symmetry. In addition, we expect that some of the flat directions may be light during 
inflation, since there are many. However, since they interact with the standard-model 
particle through the SM gauge interactions, they might be easily thermalized after the 



reheating. Interestingly, there is one possibility: Q-balls 43|]. After the flat direction starts 
oscillating, it generically experiences spatial instabilities and deforms into non-topological 
solitons called Q-balls ^ Once the Q-balls are formed, their life time is typically very 
long since the decay and evaporation processes occur only around the surfaces of the Q- 



balls 4^. The relation (|30|) can be satisfied for certain values of the inflation scale and 
the reheating temperature. For instance, we take Hi = 10^^ GeV. Assuming the gravity 



mediation or the anomaly-mediation 45[ , the typical charge of the Q-ball (of the gravity- 
mediaiton type) is Q ~ 10^^ 46|. Such Q-balls decay or evaporate depending on the 
details of the flat directions, but the typical decay or evaporation temperature is between 
C(l) GeV to C(IOO) GeV. Then the relation §^ is met for Tn ~ 10^ ~ lO^GeV. What 
makes this scenario slightly complicated is that the decay temperature is related to the 
charge of the Q-balls, which in turn depends on the initial amplitude. Therefore, in this 
model using the Q-balls, we also have to take account of the non-Gaussianity produced by 
the so-called modulated reheating scenario. Although we expect that the results derived in 
the previous section can be qualitatively valid even in this case, there must be quantitative 
difference. The detailed analysis of this scenario is beyond the scope of this paper, and 
we leave it for the future work. Nevertheless we would like to emphasize that the flat 
directions in SSM can be ungaussitons, which makes us to believe that the ungaussiton 
mechanism is indeed feasible and that one can find further explicit models. 

Each scenario described above has its own implications on the inflationary scale, the 
mass scale and the lifetime of the ungaussiton. We will discuss this issue in a separate 
paper. 



Note that the spatial seale of the Q-balls are much smaller than the cosmological scales of our interest. 
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IV. DISCUSSION AND CONCLUSIONS 



We would like to emphasize here that the ungaussiton is different from a curvaton 47 [ 



which has 



49 



50 



51 



)een extensively studied as a possible source for the non-Gaussianity [17 
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m. 



Q 



54| . The original motivation of the curvaton scenario was to explain 
the adiabatic density fluctuations by the fluctuations in the curvaton fleld. Therefore it 
is crucial for the curvaton to develop a large non-zero expectation value during inflation. 
The large expectation value helps the curvaton to dominate the energy density of the 
universe, or at least it makes the energy density of the curvaton sizable if it does not 
dominate. At the same time, the fluctuations in the energy density of the curvaton fleld 
becomes almost Gaussian with relatively small non-Gaussianity. On the other hand, the 
ungaussiton has a classically vanishing fleld value during inflation. Therefore its energy 
density exhibits strong non-Gaussianity, and it is not responsible for the almost Gaussian 
curvature perturbations observed in the CMB. Note that this difference results in the 
smaller number of the parameters in the ungaussiton scenario relative to the curvaton 
scenario, which makes the former more predictive. This is because the amplitude of the 
ungaussiton is determined solely by the Hubble parameter during inflation, while the 
amplitude of the curvaton is a free parameter in the curvaton scenario ™. Thanks to 
this feature, we have obtained the consistency relation between the bispectrum and the 
trispectrum (13T|) . Although the techniques to calculate the non-Gaussianity are common 
to the case of the curvaton, the background philosophy is different: the presence of the 
symmetry-enhanced points in the fleld space is the key idea for our arguments. 

Lastly, let us discuss how the ungaussiton can remain around the minimum of the 
potential during inflation. Since the mass is assumed to be lighter than the Hubble pa- 
rameter, one may wonder why the ungaussiton can flnd its origin and sit there during and 
after inflation. Our answer is as follows. 



In the ungaussiton scenario, there is numerical uncertainty of order unity in the effective aniphtude of 
a, since it arises from the purely quantum fluctuations. Still, we have less number of the dimensionful 
parameters. 
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There are several ways to infer how the initial position is determined, although we 
cannot decisively predict it unless the dynamics before (or during) inflation is unravelled. 
If the position is determined in a probabilistic way, the probability distribution must 
have its extremum at the origin, if the theory respects the symmetry around the origin. 
(This probability distribution has nothing to do with the quantum fluctuations. The 
initial condition we are talking here is about the position before the ungaussiton starts to 
fluctuate.) We do not know if the extremum is either (local) minimum or maximum. If 
the probability becomes maximum at the origin, the scalar field likely sits at the origin 
at the beginning of the universe. Such scalar can be an ungaussiton. On the other hand, 
the symmetry would be spontaneously broken if the probability is minimum at the origin, 
which is more suitable for the curvaton scenario. 

One can also imagine that the initial position is determined by some dynamics. For 
instance, if there is another inflationary phase before the primordial inflation, the scalar 
field may acquire a positive mass of the order of the Hubble parameter around the origin. 
Then the scalar will settle down at the origin during the pre-inflation. Or, there might be 
a radiation-dominated phase after the pre-inflation and before the last inflation starts. If 
the ungaussiton is thermalized during this phase, it will stay near the origin when the last 
inflation starts. 

If the initial position is somehow set around the origin as described above, the ungaus- 
siton will start to fluctuate about the origin during the last inflation. If it lasted for 50 — 60 
e-foldings, it still remains around the origin. However, if the last inflation lasted for much 
longer period, the situation changes. 

As the inflation lasts longer, the fluctuations will be accumulated, and the initial con- 
dition set before the inflation becomes less important, and finally, it will be forgotten in 
the end. The asymptotic distribution is known as the BD distribution. If it is reached, 
the ungaussiton takes a value of 0{Hj/m) in most regions of the whole universe. Our 
observable universe may happen to be in a (small) region where the ungaussiton takes a 
value of 0{Hi). This is a fine-tuning li m <^ Hj, since such regions are quite rare |l2l |. 
However, this is not the case if the effective mass m is not much smaller than Hj. If we 
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take m = O.lHj, for instance, the RMS of the field values is O{10Hj). Therefore, we can 
find a (sub) universe where the spatial average of the ungaussiton field is 0{Hi) without 
severe fine-tunings. Note that the fluctuations of the ungaussiton in the observable uni- 
verse are not damped significantly, for m < O.lHj, because the suppression factor is given 

2 

by ~ exp(— ^A^*) where A^* is the e-folding number from the time when the mode crossed 
the horizon to the time of the inflation end. Such possibility can be naturally realized in 
the context of supergravity. In supergravity, a scalar field generically acquires a mass of 
0{Hj) during inflation. Depending on a precise form of the Kahler potential, the effective 
mass can be lighter than Hj. (One can even make the correction vanish if the Kahler 
potential is finely tuned.) Therefore, without severe fine-tunings, it is possible that the 
effective mass is only slightly smaller than Hj. 

After all, we cannot tell the initial position of a scalar field due to our ignorance 
regarding the scalar dynamics in the inflationary phase. However, without any theoretical 
prejudice, we can still expect that a scalar field possibly sits at the origin, simply because 
the origin is so special in its field space. 

If there are the symmetry-enhanced points, we believe it plausible that the scalar fields 
stay at the points during and after inflation. If some of them are light during inflation, 
they acquire quantum fluctuations. Such light scalars are ungaussitons: their energy 
densities are necessarily highly non-Gaussian! Importantly, such non-Gaussianity must 
have been imprinted in the CMB spectrum. For the effects to be large enough to be 
actually measured by the ongoing and planned observations, the relations Eqs. fl28l) or 
( !30|) must be satisfied. Note that the criterion for a scalar field to be an ungaussiton is not 
strict at all. Any scalar field charged under some symmetries can become an ungaussiton, 
if it is light during inflation. If such scalars are mediocre in nature, we may well expect 
that there are indeed such ungaussitons satisfying the relations Eq. (!28l) or (130|) . 

We do not claim that the ungaussiton mechanism is superior to the curvaton scenario. 
These two are independent and interesting possibilities to generate sizable non-Gaussianity 
in the CMB spectrum, and as we have mentioned before, the ungaussiton scenario predicts 
distinctive bispectrum and trispectrum: the enhancement at the small scales and the 
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consistency relation fl3ip . Therefore, the future observations may be able to tell one from 
the other. 
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